PROPAGATION OF LOW REGULARITY FOR SOLUTIONS OF 
NONLINEAR PDES ON A RIEMANNIAN MANIFOLD WITH A 
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Abstract. Following [8], we introduce a notion of para-products associated to a semi-group. 
We do not use Fourier transform arguments and the background manifold is doubling, endowed 
with a sub-laplacian structure. Our main result is a paralinearization theorem in a non-euclidean 
framework, with an application to the propagation of regularity for some nonlinear PDEs. 
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The theory of paradifferential calculus was introduced by Bony in [9] and developed by many 
others, particularly Meyer in [28]. This tool that arises is quite powerful in nonlinear analysis. 
The key idea relies on Meyer's formula for a nonlinarity F(f) as M(x, D)f+R where F is smooth 
in its argument(s), / belongs to a Holder or Sobolev space, M(x,D) is a pseudodifferential 
operator (depending on /) of type (1, 1) and R is more regular than / and F(f). This operation 
is called the "paralinearization". 

Such an approach has given many important results (or improvements of existing results): 
Moser estimates, elliptic regularity estimates, Kato-Ponce inequalities, ... and is the basis of 
microlocal analysis. 

The notion of paradifferential operators is built on appropriate functional calculus and sym- 
bolic representation, available on the Euclidean space. The Fourier transform is crucial by this 
point of view to study and define the symbolic classes. That is why this approach cannot be 
extended to Riemannian manifolds. 
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However, for the last years, numerous works deal with nonlinear PDEs on manifolds. So it 
seems important to try to extend this tool of "paralinearization" in a non-Euclidean situation. 

First, on specific situations, namely on a Carnot group it is possible to define a suitable Fourier 
transform, involving irreductible representations. In this context, we can also define the notion 
of symbols and so of pseudo-differential calculus (see the survey [5] of Bahouri, Fermanian- 
Kammerer and Gallagher for Heisenberg groups and [21] of Gallagher and Sire for more general 
Carnot groups). Excepted this particular setting, no Fourier transform are known. 

Following this observation, the aim of this current work is to define another suitable notion 
of paralinearization on a manifold, without requiring use of Fourier transform. Since (non- 
linear) PDEs on a manifold usually requires vector fields, we work on a manifold having a 
sub-Riemannian structure. To define a suitable paralinearization, we use paraproducts defined 
via the heat semigroup (introduced by Bernicot in [8], independently by Frey in [19, 20] and 
already used by Badr, Bernicot and Russ in [4] to get Leibniz type estimates and algebra prop- 
erties for Sobolev spaces) and look for a paralinearization result. However, a new phenomenom 
appears (due to the lack of flexibility of the method), the classical paralinearization result holds 
only for low regularity. More precisely, we prove the following (see Section 3 for Sobolev spaces 
and Section 4 for the definition of the paraproduct n): let consider Sobolev spaces associated to 
a Sub-Laplacian operator on a Riemannian manifold then (under usual assumptions), we have 

Theorem 0.1. Consider p G (l,oo), s G (d/p,l) and f G W s+e ' p for some e > (as small as 
we want). Then for every smooth function F G C°°(R) with F(0) = 0, 

(1) F(f)=U F , (f) (f) + w 

with w g w 2s ~ d / p ' p . 

As in the Euclidean situation, we are able to obtain some applications concerning propagation 
of the regularity for solutions of nonlinear PDEs. 

With respect to the well-known paralinearization results, the first point is that we have only 
a gain of regularity at order s — d/p — e and the main difference is that this result is only proved 
for s < 1. This condition can be viewed as very strong but we will explain in Remark 5.4 how to 
allow larger regularity s > 1 with adding extra terms in the paraproducts (involving the higher 
order derivatives of the nonlinearity F). This limitation s < 1 already appeared in [14] and [4]: 
the use of higher order Sobolev spaces require to understand higher order Riesz transforms and 
some cancellations properties of the iterated laplacians A k ... which seems to be very difficult. 



1. Preliminaries : Riemannian structure with a sub-Laplacian operator 

In this section, we aim to describe the framework and the required assumptions, we will use 
after. Let us precise the main hypothesis about the manifold M and the operator L. 

1.1. Structure of doubling Riemannian manifold. In all this paper, M denotes a complete 
Riemannian manifold. We write \i for the Riemannian measure on M, V for the Riemannian 
gradient, | • | for the length on the tangent space (forgetting the subscript x for simplicity) and 
|| • \\lp for the norm on L p := L p (M,fi), 1 < p < +oo. We denote by B{x,r) the open ball of 
center x G M and radius r > 0. 

1.1.1. The doubling property. 

Definition 1.1 (Doubling property). Let M be a Riemannian manifold. One says that M 
satisfies the doubling property (D) if there exists a constant C > 0, such that for allx G M, r > 
we have 



fj,(B(x,2r)) < Cfi(B(x,r)). 
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Lemma 1.2. Let M be a Riemannian manifold satisfying (D) and let d := Zo^C. Then for all 
x, y 6 M and 9 > 1 

(2) K B ( X , eR )) < Ce d n(B{x, R)) 
There also exists c and N > 0, so that for all x,y € M and r > 

(3) »(B(y, r))<c(l + " ^(B(x, r)). 



For example, if M is the Euclidean space M = M. d then N = and c = 1. 
Observe that if M satisfies (D) then 

diam(M) < oo 44> /i(M) < oo (see [1]). 

Therefore if M is a complete Riemannian manifold satisfying (D) then fJ.(M) = oo. 

Theorem 1.3 (Maximal theorem). (|ll]j Let M be a Riemannian manifold satisfying (D). 
Denote by A4 the uncentered Hardy-Littlewood maximal function over open balls of M defined 
by 

Mf(x) := sup \f\ Q 



Q ball 



where fi 



:= + fdfi := f fdf-L. Then for every p € (l,oo], M. is L p bounded and moreover 

Je V( e ) Je 

of weak type (1,1). 

Consequently for s £ (0, oo), the operator A4 S defined by 

M s f(x) := [M(\f\ s )(x)} 1/S 
is of weak type (s, s) and LP bounded for all p € (s, oo]. 

Doubling property allows us to control the growth of ball-volumes. However, it can be inter- 
esting to have a lower-bound too. So we will make this following assumption: 

Assumption 1.4. We assume that there exists a constant c > such that for all x € M 

(4) fi(B(x, 1)) > c. 

Due to the homogeneous type of the manifold M , this is equivalent to a below control of the 
volume (MVd) 

(MVd) fi(B(x,r))>r d 

for allO <r < 1. 

1.1.2. Poincare inequality. 

Definition 1.5 (Poincare inequality on M). We say that a complete Riemannian manifold M 
admits a Poincare inequality (P q ) for some q G [1, oo) if there exists a constant C > such 
that, for every function f € W^(M) (the set of compactly supported Lipschitz functions on M) 
and every ball Q of M of radius r > 0, we have 

(P q ) (J 1/ - fo\ q d^' q <Cr (J iV/lW) 1 9 . 

Remark 1.6. By density ofC^(M) in W^(M), we can replace W^(M) by C7 °°(M). 

Let us recall some known facts about Poincare inequalities with varying q. 
It is known that (P q ) implies (P p ) when p > q (see [24]). Thus, if the set of q such that (P q ) 
holds is not empty, then it is an interval unbounded on the right. A recent result of S. Keith 
and X. Zhong (see [26]) asserts that this interval is open in [1, +oo[ : 
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Theorem 1.7. Let (X,d,fi) be a complete metric-measure space with fi doubling and admitting 
a Poincare inequality (P q ), for some 1 < q < oo. Then there exists e > such that (X,d,/x) 
admits (P p ) for every p > q — e. 

Assumption 1.8. We assume that the considered manifold satisfies a Poincare inequality (Pi). 
Indeed we could just assume a Poincare inequality (Pa) for some a < 2 and all of our results 
will remain true for Lebesgue exponents bigger than a 

1.2. Framework of Semigroup. Let us recall the framework of [17, 18]. 
Let cj G [0, 7r/2). We define the closed sector in the complex plane C by 

S u := {zeC, |arg(*)| <^}U{0} 

and denote the interior of S u by 5°. We set i? 0O (S'£) for the set of bounded holomorphic 
functions b on S®, equipped with the norm 

WHh^so) : = IHU°°(s°)- 

Then consider a linear operator L. It is said of type cj if its spectrum o~(L) C S w and for each 
v > oj, there exists a constant c„ such that 

Hex- a)- 1 !!^^ < c^iai- 1 

for all A ^ S v . 

We refer the reader to [17] and [27] for more details concerning holomorphic calculus of such 
operators. In particular, it is well-known that L generates a holomorphic semigroup (A z := 
e _zL ) zg 5 jr/2 _ u . Let us detail now some assumptions, we made on the semigroup. 

Assume the following conditions: there exists a positive real m > and 5 > 1 with 

Assumption 1.9. • For every z € S n /2-u, the linear operator A z := e~ zL is given by a 

kernel a z satisfying 

1 / d(x vV - d - 2N ~ s 

(5) Mx,y)\< , n/ ^ (l + 



n(B(x,\z\ 1/2 )) \ \z\ 1/2 
where d is the homogeneous dimension of the space (see (2)) and N is the other dimension 
parameter (see (3)); N > could be equal to 0. 

• The operator L has a bounded H ^-calculus on I? . That is, there exists c u such that for 
b € iT 00 (5°) ; we can define b(L) as a I? -bounded linear operator and 

(6) \\b(L)\\ L 2^ L 2 < evII&Hoo. 

• The Riesz transform 1Z := VL" 1 / 2 is bounded on LP for every p S (1, oo). 

Remark 1.10. The assumed bounded Hoc-calculus on L 2 allows us to deduce some extra prop- 
erties (see [18] and [27]) : 

• Due to the Cauchy formula for complex differentiation, pointwise estimate (5) still holds 
for the kernel of (tL) k e~ tL with t > 0. 

• For any holomorphic function ip S H(S®) such that for some s > 0, ^5 i+|zp s ' ^ e 
quadratic functional 

/ roc 1/2 

/-►( / imm 2 



IS 



L? -bounded. 



Remark 1.11. Concerning a square estimate on the gradient of the semigroup, it follows that 
for every integer k > the square functional 

I \t^ 2 V(tL)^ L (f)\ d ± 
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is bounded on I? . Indeed, this is just a direct consequence of the boundedness of the previous 
square functions, by making appear the Riesz transform and uses its L 2 -boundedness. 

Remark 1.12. We claim that Assumption (7) is satisfied under the L 2 -boundedness of the Riesz 
transform 1Z : = VL -1 / 2 . 

Indeed iflZ is L 2 -bounded, then it admits L 2 -valued estimates, which yield 



2 



dtdfi\ 1/2 ^ ( f°° fe+1/2 tL 2 dtdfi\ 



This gives the desired result 



1/2 



^| ( '/ M <L)'e-<V)| 2 ^)' /! 



< 



l/l 



L- ! 



thanks to Remark 1.10. 

About different square functions, we have the following proposition: 

Proposition 1.13. Under these assumptions, we know that the square functionals in Remark 
1.10 or in (7) are L p -bounded for every p € (l,oo). 

Proof. Let T be one of the square functions in Remark 1.10. We also already know that it is L 2 
bounded, by holomorphic functional calculus. Then consider the "oscillation operator" at the 
scale t: 

Bt := 1 - At = 1 - e~ tL = - [ Le~ sL ds. 



Then, by using differentiation of the semigroup, it is classical that TBt satisfies L 2 — 1? off- 
diagonal decay at the scale i 1//?n , since the semigroup e~~ tL is bounded by Hardy Littlewood 
maximal function A4 S _. So we can apply interpolation theory (see [7] for a very general exposi- 
tion of such arguments) and prove that T is bounded on IP for every p € (s_, 2] (and then for 
p € [2, oo) by applying a similar reasoning with the dual operators). 

Then consider a square function U of type (7). Then by using the Riesz transform, it yields 

/ roo ja 1/2 

v(f) = \J \mtL)f\ 2 - 

with ip(z) = z l / m 4>{z). Since 1Z is supposed to be L p -bounded, it verifies £ 2 -valued inequalities 
and so the L p -boundedness of U is reduced to the one of a square functional of previous type, 
which was before proved. □ 

1.3. Framework of sub-Laplacien. We will only consider operators L which are sub-Laplaciens 
and generating a semigroup (e~* L )t>o satisfying the above assumptions. Let us first explain what 
a sub-Laplacian means : 

We assume that there exists X = {Xk}k=i,...,K a finite family of real-valued vector fields (so 
Xk is defined on M and Xk{x) G TM X ) such that 

(8) £ = -XX 

k=l 

We identify the X^s with the first order differential operators acting on Lipschitz functions 
defined on M by the formula 

X k f(x) = X k (x) ■ V/(x), 
and we set Xf = (X ± f, X 2 f, • • • , X K f) and 

\Xf(x)\=(jp\X k f(x)\* S j , x E M. 
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We define also the higher-order differential operators as follows : for / C {1, n} k , we set 

iel 

We assume the following: 

Assumption 1.14. For every subset I , the Ith-local Riesz transform TZj := Xj(l + L)~^l/ 2 and 
its adjoint TV\ := (1 + L)~I 7 I/ 2 X/ are bounded on LP for every p G (1, oo). 

Remark 1.15. It is easy to check that this last assumption is implied by the boundedness of 
each local-Riesz transform IZi and 1Z* in Sobolev spaces W k,p for every p G (1, oo) and k G N. 
Indeed for I = {i±, ...,i n }, we have 

m/WLP < \\Ki(x i2 ...x in f)\\ w m-^ z mL...,ijh P . 

Repeating this reasoning, we obtain that the Sobolev boundedness of the Riesz transforms and of 
its adjoint implies the previous Assumption. 

From now on, we will consider a doubling Riemannian manifold M satisfying Poincare in- 
equality (-Pi), lower bound of the volume Assumption (1.4) and a structure of sub- Riemannian 
laplacian associated to a semigroup satisfying Assumption (1.9) and with bounded Riesz trans- 
forms (Assumption (1.14)). 

2. Examples of such situations 

In this section, we would like to give two examples of situations where all these assumptions 
are satisfied. 

2.1. Laplacian operators on Carnot-Caratheodory spaces. Let Q be an open connected 
subset of M. d and Y = {Yjfc}fc =1 a family of real- valued, infinitely differentiable vector fields. 

Definition 2.1. Let f2 and Y be as above. Y is said to satisfy Hormander's condition in f2 if 
the family of commutators of vector fields in Y (Yi, [Yi,Yj], ....) span M. d at every point of ft. 

Suppose that Y = {Yk}^ =1 satisfies Hormander's condition in Q. Let Cy be the family of 
absolutely continuous curves C : [a, b] — > £1, a < b, such that there exist measurable functions 
Cj (t), a < t < b, j = l,--- ,M, satisfying ZjLi^t) 2 < 1 and ?(t) = Zf=i <H (<)^(C(<)) for 
almost every t G [a, b]. If x, y 6 define 

p(x,y) = inf{T > : there exists £ G Cy with ("(0) = x and COO = !}• 

The function p is in fact a metric in called the Carnot-Caratheodory metric associated to Y. 
This allows us to equipped the space of a sub- Riemannian structure. 

2.2. Lie groups. Let M = G be a unimodular connected Lie group endowed with its Haar 
measure d\i = dx and assume that it has a polynomial volume growth. Recall that "unimodular" 
means that dx is both left-invariant and right-invariant. Denote by C the Lie algebra of G. 
Consider a family X = {X X K } of left-invariant vector fields on G satisfying the Hormander 
condition, which means that the Lie algebra generated by the AYs is C By "left-invariant," 
one means that, for any g G G and any / G Cq°(G), X(r g f) = r g (Xf), where r g is the left- 
translation operator. As previously, we can build the Carnot-Caratheodory metric on G. The 
left invariance of the Xj's implies the left-invariance of the distance d. So that for every r, 
the volume of the ball B(x,r) does not depend on x G G and also will be denoted V(r). It is 
well-known (see [23, 29]) that (G,d) is then a space of homogeneous type. Particular case are 
Carnot groups, where the vector fields are given by a Jacobian basis of its Lie algebra and satisfy 
Hormander condition. In this situation, two cases may occur : either the manifold is doubling 
or the volume of the balls admit an exponential growth [23]. For example, nilpotents Lie groups 
satisfies the doubling property ([16]). 
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We refer the reader to [30, Thm5.14] and [14, Section 3, Appendix 1] where properties of 
the heat semigroup are studied: in particular the heat semigroup e~ tL satisfies gaussian upper- 
bounds and Assumption 1.9 on the higher-order Riesz transforms (Assumption (1.14) is satisfied 
too. 

2.3. Carnot groups. Particular cases of nilpotents Lie groups are the Carnot groups (if it 
admits a stratification). A stratification on a Lie group G (whose q is its Lie algebra) is a 
collection of linear subspaces V±, V r of G such that 

g = Vi © ... © V r 

which satisfy [Vi, Vi] = Vf+i for i = 1, r — 1 and [V\, V r ] = 0. By [Vi, Vi] we mean the subspace 
of G generated by the elements [X, Y] where X G Vi and Y G V^. Consider rtj the dimension of 
Vi, n\ + • • • + n r = d and dilations {5a}a>o of the form 

S x ( x ) = (AxW,AV 2) ,-- - A s x {r) ), G V. 

The couple G = (G, 5\) is said to be a homogeneous Carnot group (of step r and n\ generators) 
if S\ is an automorphism of G for every A > and if the first n\ elements of the Jacobian basis 
of g, say Zi, ■ ■ ■ , Z m , satisfy 

(9) rank(Lie[Zi, • • • , Z ni ](x)) = d, for all x G G, 

where Lie[Zi, • • • , Z ni ] is the Lie algebra generated by the vector fields Z±, ■ ■ ■ Z ni . The number 
Q = J2i=i i n i i s called the homogeneous dimension of G. 

As for example the different Heisenberg groups, T~i d is a Carnot group of dimension Q = 2d+2. 
We refer the reader to [21] for an introduction of pseudo differential operators in this context 
using a kind of Fourier transforms involving irreductible representations (and to [5] for a complete 
work about pseudo-differential calculus on Heisenberg groups). 

2.4. Riemannian manifolds with a bounded geometry. We shall say that a Riemannian 
manifold M has a bounded geometry if 

• the curvature tensor and all its derivatives are bounded 

• Ricci curvature is bounded from below 

• and M has a positive injectivity radius. 

In such situations, we know that there exists a collection of smooth vector fields Xi,...,X K 
such that 

i=l 

Moreover Assumptions 1.9 and 1.14 are satisfied (see [14] and [34]). 

3. The scale of Sobolev spaces 
We use the Bessel-type Sobolev spaces, adapted to the operator L: 
Definition 3.1. For p G (1, oo) and s > 0, we set 

W s, P = w w := l f e LP) (! + L) s ' 2 {f) G L p } . 

First, we have this characterization: 
Proposition 3.2. For all p G (1, oo) and s > 0, we have the following equivalence 

\\f\\ LP + \\L^(f)\\ LP ^\\(i + Ly/ 2 f\\ LP . 
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Proof. Set a = s/2 and write a = k + 9 with k £ N and 9 G [0, 1). We decompose (1 + L) a with 
the semigroup as following 

f°° dt 

(1 + L) a f = / e-'e-* L (l + L)t l - 6 ^(l + L) fc (/) 
J o 1 

= e - t e- tL t 1 - e -(l + L) k (f)+ H e- t e-*(tL) x -°-L 9 (l + L) k {f). 
Jo t Jo * 

The first integral-operator is easily bounded on LP since the semigroup e~ tL is uniformly 
bounded. The second integral operator is bounded using duality: 

( /V^W^WpS) = r \^(e^'\tL)¥(u),e- a ''\tL*)¥g)*± 



o 



< 



DC 



1/2 



oo 



2 dt 



1/2 



Since (1 — a)/2 > 0, then the two square functionals are bounded in LP and LP (by Proposition 
1.13) and that concludes the proof of 

11(1 + LYf\\ LP < ||(1 + L) k f\\ LP + \\L 9 (1 + L) k (f)\\ LP . 

Then, developing (1 + L) k , it follows a finite sum of ||L z (/)||j> with z € [0,a]. We decompose 



f°° rlt f' 1 rlt f' 00 rlt 

/ e- tL {tL)°U)t-'^ r = e - tL {tLYU)r^+ e ~ tL (tL) a t^. 
Jo * jo 1 Ji 1 



The first quantity in LP is controlled by ||L Q (/)||i P and the second one by ||/||lp, which concludes 
the proof of 

ll(l + ^r /2 /l|L P <||/||L P + ||^ /2 (/)||L P . 

Let us now check the reverse inequality. As previously, for u = or u = a we write 

poo f lf 

L u f= / e - t{l+L) (l + L)L u t 1+a — (l + L) a f. 
Jo t 

By producing similar arguments as above, the operator J °° e - *( 1+L )(l + L)L u t 1+OL Y is easily 
bounded on LP (splitting the integral for t < 1 and t > 1) and we can also conclude to 

\\L u (f)\\ LP <\\(i+Lrf\\ LP , 

which ends the proof. □ 
Corollary 3.3. For all p € (1, oo) and < t < s, we have the following inequality 

||^/IIlp<II(i + ^)VI|lp^||/|| w **. 

Let us then describe classical Sobolev embeddings in this setting (see [4] for a more general 
framework) : 

Proposition 3.4. Under Assumption 1.4 (lower bound on the ball-volumes), let s > t > be 
fixed and take p < q such that 

1 t 1 s 

q d p d 
Then, we have the continuous embedding 

We refer the reader to [4, Proposition 3.3] for a precise proof. The proof is based on a 
spectral decomposition, to write the resolvant with the semigroup and then to use the off- 
diagonal estimates (here the pointwise estimates on the heat kernel). 

Corollary 3.5. Under the previous assumption, W s,p L°° as soon as 

d 

s > -. 
p 



PARA-DIFFERENTIAL CALCULUS ON A MANIFOLD 9 

We now recall a result of [4], where a characterization of Sobolev spaces is obtained, involving 
some fractional functionals. 

Proposition 3.6. [4, Thm 5.2] Under Poincare inequality (Pi), for s G (0,1) we have the 
following characterization : a function f G LP belongs to W s ' p if and only if 

SP/ oc f(x) 



i / r Y /p 




r s yfi(B(x,r)) J B ( Xir ) 

belong to LP , for some p < min(2,p). 

This characterization can be extended for s > 1, using the sub-Laplacian structure. Indeed, 
we have this first Lemma: 

Lemma 3.7. For every integer k and p G (1, oo), 

\\f\\ W H, P C* Yl W X l(f)\\Lr. 

IC{l,...,K} k 

Proof. As point out in [14], this is consequence of Assumption 1.14 about the local Riesz trans- 
forms. Indeed, for k > 1 and L a subset, we have assumed that the 7-th Riesz transform TZj are 
bounded on L p , which is equivalent to 

ll*/(/)Hi> < H/llwrl/IAp. 

Moreover, making appear the adjoint of the Riesz transforms and the resolvant (which are all 
bounded on L p ) as follows 

ft K 

(1 + L) 1 / 2 = (1 + L)-V2(i + L ) = (l + L)-^ 2 + ^2(l + L)-^ 2 Xf = (l + LyW + ^iqXi, 

i=l i=i 
we conclude to the reverse inequality and so we have proved the desired result for k = 1. We 
let the details for k > 2 to the reader, the reasoning is exactly the same but technically to be 
written (we make appear a finite sum of higher order Riesz transforms ...) □ 

We also deduce the following characterization (see Proposition 19 in [14]): 

Proposition 3.8. Let s := k + t > 1 (with k an integer and t G (0, 1) ), then 

f G W s ' p ^feL p and VI C {1, K } k , X/(/) G W l ' p 

^f£L p and VI C {1, K } k , Sf (*/(/)) G LP. 

We also deduce the following chain rule (see Theorem 22 in [14] for a proof by induction on 

k): 

Proposition 3.9. // F G C°° with F(0) = and let s := k + t > ^ (with p G (l,oo), k an 
integer and t G (0, 1) ). Then 

ITO)lk^<ll/lk- + II/IIW 

Lf F(0) ^ 0, we still have such inequalities with localized Sobolev spaces. 

We refer the reader to [14] for a proof by induction on k. Here for completeness, we produce 
another direct proof. 

Proof. We use the previous characterization of the Sobolev space with Sf for s = k+t. First using 
the differentiation rule, it comes Xi(F(f)) = Xi(f)F'(f), then XjXi(F(f)) = XjXi(f)F'(f) + 
Xi(f)Xj(f)F"(f) ... By iterating the reasoning, for / C {l,...,n} k , estimating Xj(F(f)) in 
\yt,p j g rec iuced to estimate quantities as 

i 



h :-- 



.£*=! 



(f)F {n \f) 



10 FREDERIC BERNICOT AND YANNICK SIRE 

where i a C I , n < k and ^ \i a \ = \I\ < k. Then for x, y, we have 

Hx)-h(y)\<^\X l0 (f)(x)-X lp (f)(y)\H snp |X 1q (/)(z)|||fW(/)||^ 

P a^ Z=X ' y 

+ [] snp |X ia (/)(^)|| J pW(/)(x)-^)(/)(y)|. 

z=x,y 

a 

By this way, since p < p let us choose exponents /9 a , p a such that 

1 x 1 

~ = ' Pa<Pa 

P ^ Pa 

and 

Moreover we require that 

(10) J__M+1>I * 

p a a pa 

This is possible since ^ Q \i a \ = \I | < s — i and s > d/p (indeed we let the reader to check that 
Pa = FTC P * s a g°°d choice). Moreover, we chose exponents p^, p, p~a~ and p such that 

1 v-^ 1 1 

~=V = + =, Pa<Pa 
P « Pa P 

and p < /5 < p with 

1 1 1 

- = > — + -. 

p ^ p a p 

As previously, we require (10) with p^ instead of p a and 

lis 

11 ->--!■ 

p p a 

Such exponents can be chosen by perturbing the previous construction with a small parameter 
since s > d/p. By this way (with Holder inequality), we deduce that 

S P t(h) II M Pa [X la (f)}\\F in \f)\\ L ™ 

P a+p 

+ \{M^[X la {f)]S?{F^{f). 

a 

Since F is supposed to be bounded in C°°, then F^ is Lipschitz and so, we finally obtain 
S?(h) J] M Pa [X la (f)}\\F [n \f)\\L- 

P a^p 

+ \{M^[X la if)]S p t if). 

a 

Then applying Holder inequality, we get 

\mh)\w £EIK"(^(/))IU II \\M Pa \XiM)]\\ L p a ii^ {n) (/)ii^ 

P a^P 

+nii- M ^ x ^(/)]iiL^||w)|| LF - 

a 

Since (10) with Sobolev embeddings (Proposition 3.4), we have 

\\M Pa [X ta (f)\\\ LPa < \\f\\wM*~ Sll/llw* 

and 

\\ S t P ( X ip(f))\\ L Pp ~ Wf\\ w \*i3\+t,Pfj < \\f\\w°<p- 
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So with (10) and (11), we finally obtain 

\\Snh)\\LP<\\f\\w-* + \\f\\Uv, 

where we used s > d/p and the Sobolev embedding W s,p C L°° with the smoothness of F to 
control ||F( n )(/)|| L oo. 

Since F(0) = and F is Lipschitz, we also deduce that F(f) belongs to IP, which allows us 
to get the expected result 

\\F(f)\\w^ < \\F(f)\\ LP + £ \\S?(F(f))\\ LP < \\f\\ w .„ + \\f\\ k Ws , P . 

i 

□ 

Remark 3.10. If F G C°° with F(Q) = and s > d/p then to obtain 

\\Hf)\\w-*<\\f\\w-* + \\f\\$r.*, 

it is sufficient to assume that F is locally bounded in C°° and then the implicit constant will 
depend on ||/||z,°°. Indeed, using Sobolev embedding, we know that as soon as s > d/p, W s ' p is 
continuously embedded in L°° . 

4. Paraproducts associated to a semigroup 

Our aim is to describe a kind of "paralinearization" results. In the Euclidean case, this is 
performed by using paraproducts (defined with the help of Fourier transform). Here, we cannot 
use such powerful tools, so we require other kind of paraproducts, defined in terms of semigroup. 
These ones were introduced by the first author in [8], already used in [4] and more recently was 
extended in [19, 20]. Let us recall these definitions. 

4.1. Definitions and spectral decomposition. We consider a sub-Laplacian operator L sat- 
isfying the assumptions of the previous sections. We write for convenience cq for a suitably 
chosen constant, ip(x) = cqx n e~ x {1 — e~ x ) and so 

ML) :=c (tL) N e- tL (l-e- tL ), 
with a large enough integer N > d/2. Let <j) be the function 

roc 

<f>(x) := -co / y N e-y(l - e~y)dy, 
£(*):= -co/ y N ^e^(l-e-y)dy, 



and set 4>t{L) := <f)(tL). Then we get a "spectral" decomposition of the identity as follows 
(choosing the appropriate constant cq), we have 

,dt 

T' 



r°° t 

f = - <//{tL)f- 
So for two smooth functions, we have 



fg := - [ <t>'(sL) [(/>'(uL)f <f>'(vL)g] 

Js,u,v>0 suv 

Since 4>'(x) = if)(x) := cqx n e~ x (l — e~ x ) and xcf>'{x) = 4>'(x), it comes that (by integrating 
according to t := min{s, u, v}) 

dt 



fg:=- I tP(tL) <j>{tL)f<l>{tL)g f - / 4>{tL) ${tL)f cj>{tL)g 



(12) 



<t>(tL) cj){tL)fij{tL)g 



t 
dt 
T' 



t 
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Let us now focus on the first term in (12) : 



Hf,g) 



i>(tL) <j>(tL)f<l>(tL)g 



dt 

T' 



Since N >> 1, let us write ip(z) = zip with ip (still vanishing at and at infinity). Then using 
the structure of the sub-Laplacian L, the following algebra rule holds 

L(fg)=L(f)g + fL(g) + (Xf-Xg), 

where X is the collection of vector fields Xf := (X\f, ...,X K f). Hence, we get 

dt 



i(f,g) 



+ 



+ 



4(tL)(tL) [<t>{tL)f <j){tL)g 
$(tL) \tL4>{tL)f4>{tL)g 

3 r~ 

$(tL) <t>(tL)ftI4(tL)g 



t 



dt 

T 

dt 

T 



$(tL)t(X0(tL)f ■ X0(tL)g)j. 



Combining with (12), we define the paraproduct as follows : 

Definition 4.1. With the previous notations, we define the paraproduct of f by g, by 

dt 



n fl (/) := 



i>{tL) tL<j>{tL)f<j>{tL)g 



t 



cp{tL) i>(tL)f <f>(tL)g 



dt 

T' 



Remark 4.2. We first want to point out the difference with the initial definition in [8]. There, 
general semigroup was considered and the previous operation on the term I can be performed by 
making appear the "carre du champ" introduced by Bakry and Emery (see [6] for details) 

F(f,g) :=L(fg)-L(f)g-fL(g) 

instead of the vector field X. However in [8], the paraproducts was only defined by the second 
term. This new definition comes from the following observation: considering the quantity /(/, g) 
and distributing the Laplacian as we have done (or make appearing the "carre du champ"), 



it comes three terms. The term ip(tL) tL(p(tL)f 4>{tL)g 



has the same regularity properties 



as <t>(tLi) ip(tL)f (p(tL)g (in the sense that tL<p{tL) can be considered as ip(tL)). This also 

legitimate to add this extra term in the definition of the paraproducts. 

By this way, as we will see in the next properties and in Remark 5.4, this new paraproduct is 
the "maximal" (in a certain sense) part of the product fg, where the regularity is given by the 
regularity of f . 

It naturally comes the following decomposition : 

Corollary 4.3. Let f, g be two smooth functions, then we have 

fg = U g (f) + U f (g) + Rest(f,g) 

where the "rest" is given by 



Rest(f,g) :-- 



^tL)(t^ 2 X^tL)f, t l l 2 X<jy{tL)g) 



dt 
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4.2. Boundedness of paraproducts in Sobolev and Lebesgue spaces. Concerning esti- 
mates on these paraproducts in Lebesgue spaces, we refer to [8]: 

Theorem 4.4 (Boundedness in Lebesgue spaces). Forp,q G (l,oo] with < £ := | + | then 

(f,g)^U g (f) 

is bounded from LP x L q into U . 

Let us now describe boundedness in the scale of Sobolev spaces. 

Theorem 4.5 (Boundedness in Sobolev spaces). For p,q,r G (l,oo) with - := - H — and 
s G (0; 27V - 4) then 

(f,g)^U g (f) 

is bounded from W s;p x L q into W s ' r . 

Proof. It is sufficient to prove the following homogeneous estimates : for every ft G [0, N — 2) 

||^(/)||^<||^(/)IMMIl«. 

For f3 = 0, this is the previous theorem so it remains us to check it for ft G (0, TV — 2). We recall 
that 



n 9 (/) 



^(tL) tL<j>(tL)f<f>(tL)g 



dt 

T 



<K*L) t/>(tL)f <f>(tL)g 



dt 



giving rise to two quantities, n*(/) and H 2 g {f). Indeed, applying L 13 to the paraproduct Il 2 (f), 
it yields 



L^(f) 



L?4>{tL) [^(tL)f<P(tL)g 

l 

$(tL) [t-fy(tL)f<l>(tL)g 
^(tL) ^{tL)L? f <j>(tL)g 



dt 

T 

dt 

T 

dt 

T' 



where we set VK- 2 -) = z l3 (f)(z) and ^(z) = z ^^(z). So if the integer iV in and ^ is taken 
sufficiently large, then ^ and ip are still holomorphic functions with vanishing properties at 
and at infinity. As a consequence, we get 

L^U g (f)=U 9 (L^f) 

with the new paraproduct II built with ip and ip. We also apply the classical reasoning aiming 
to estimate this paraproduct. By duality, for any smooth function h G L r we have 



<^n? (/),/»> 



^(tL*)hi>(tL)(LPf)j(tL)gjd» 



< 



mtL*)h\ 2 ^ 



■if 

From the pointwise decay on the semigroup (5), we know that 



\i/>(tL)(lff)\ 2 ^ ) sup|0(tL) 5 |d/i. 



sup|0(tL) 5 (x)| <M(^)(s) 
t 



and so by Holder inequality 



mL*)h\ 2 



dt 



1/2 



mL)(L p f)\ 



, dt 



1/2 



IIMsl 



LP 



14 



FREDERIC BERNICOT AND YANNICK SIRE 



Since ip and ip are holomorphic functions vanish at and having fast decays at infinity, we know 
from (1.13) that the two square functions are bounded on Lebesgue spaces. We also conclude 
the proof by duality, since it comes 



(L^U 2 Jf),h) 



< 



Iff 



LP 



\9\\li 



We let the reader to check that the same arguments still holds for the first part n*(/) and so 
the proof is also finished. □ 

5. Linearization Theorem 

Theorem 5.1. Consider s € (d/p, 1) and f € W s+e ' p for some e > 0. Then for every smooth 
function F G C°°(R) with F(0) = 0, 

(13) F(f) = U FV) (f)+w 

with w E W 2s - d / p ' p . 

We follow the proof in [10, 28, 9]. 

Proof. Let us refer the reader to the operators <f>{tL) and ip(tL), defined in Subsection 4.1: 
ijj(x) = cqx n e~ x (1 — e~ x ), (j) is its primitive vanishing at infinity. Let us write ip(z) = z~ l ijj(z) 
and (j) its primitive vanishing at infinity. Moreover, these functions are normalized by the suitable 
constant cq such that 0(0) = 1. 



It comes 



and so we decompose 



/ = Jim0(tL)(/) 



F(f) = <p(L)F^(L)f) 



1 ±MtL)FWL)f)dt. 
dt 



Since 



t-0(tL)F(<p(tL)f) = tL4>'{tL)FQ>{tL)f) + <f>(tL) [(tL#(tL)f)F'$(tL)f) 
= cj>'{tL)F${tL)f) + 4>{tL) \(<p'(tL)f)F'(0(tL)f) 



we get 



F(f) = ct>(L)F(cf>(L)f) 
= 4>{L)FQ>{L)f) 



1 ~ ~ r i dt 

(j) (tL)tL[F((f)(tL)f)\ + 4>(tL) [{$ \tL)f)F'{<t>{tL)f)) j 

1 4>'{tL) \ F" ((f)(tL) f)^ 1 / 2 X(f>(tL) f\ 2 + F'Q>{tL)f)tL4(tL)f 



+ <P(tL) ( ( /> , (tL)f)F , (<l>(tL)f 



dt 

T' 



where we used the differentiation rule for the composition with the vector fields X = (Xi, X K ) 
We also set 

w-I + II + III + IV + V 

with 



II :-- 



I := d>(L)F^(L)f), 
\tL) \F"{4>{tL)f)\ty 2 X4>{tL)ff 



o 



dt 

T' 



III :-- 

IV :-- 



1 4/{tL) [fatL)F'(f) - F'{cp{tL)f)) tL</>(tL)f] j: 

dt 

7' 



ct>{tL) (</>'{tL)f) U(tL)F'(f) - F>(0(tL)f) 



PARA-DIFFERENTIAL CALCULUS ON A MANIFOLD 



15 



and 



V :-- 



4>(tL) tL4>(tL)f4>{tL)F'(f) 



dt 



+ 



0(tL) *P(tL)f<P(tL)F'(f) 



dt 



in order that (13) is satisfied. It remains us to check that each term belongs to W 2s d l p > p . 
Step 1: Term I. 

Since / £ W s+e ' p then 4>{L)f belongs to W p ' p for every p > s + e and so Proposition 3.9 yields 
that 



<t>(L)F{<j>(L)f) 



< 



W2a-d/p,p 



Step 2: Term V. 

We only treat the first term in V (the second one can be similarly estimated). Using duality, 
we have with some g € L p ' and for a £ {0, 2s — d/p} since a > and t > 1 



< 



2 dt 



1/2 



/roc ____ 
y (tL) a / 2 ^(tL*)g tL4>(tL)f4>(tL)F\f) 

tL0(tL)fj(tL)F'(f) 
\ LP S upU(tL)F'(f)\\ L ~, 



< 



where we used the boundedness of the square functional. Then we conclude since (p(tL)F' ( f) 
is uniformly bounded by H^'C/OIl 00 which is controlled by ||/||w, P (due to Sobolev embedding 
with s > d/p and Proposition 3.9). 

Indeed our problem is to gain some extra regularity (from s to 2s — d/p) so the main difficulty 
relies on the study of the "high frequencies" and not on the lower ones. 

Step 3: Term IF 

By duality and previous arguments, we get 



\II\ 



W2s-d/p,p 



< 



Of t ~ 2s+d/p \ ti/2x ^ tL )f\ i j) 









r 


+ 


a 




LP 





dt\ 



1/2 



LP 



where we decomposed the norm with its homogeneous and its inhomogeneous parts and then 
used uniform boundedness of F" {4>{tL) f) . Since (using Z^-boundedness of the Riesz transforms, 
see Assumption 1.14 and Sobolev embedding) 

\\X4(tL)f\\\ L «> < \\X$(tL)f\\ wd/p+e , p 

< \\L^(tL)f\\ L p + \\L l ' 2+d / 2p+ ^{tL)f\\ L p 

<f /a - 1/a ||/lk.*+f /2 - i /*- 6 / a - 1 / a ||/|| W r., 
< t s/2- d /2 P -e/2-l/2 mwsp 

where we used t < 1. Finally it comes, 



HI 



\iy2s-d/p,p < 



< 



j\-"\^Xm)f\^' 2 < (jf 1 |(fL)V^( S ^)/2^ L)L ( S+ e)/2 /| 2^ 



1/2 



LP 



where we used s < 1 and the boundedness of the square functional. 
Step 4: Terms III and IV. 

For these terms, we follow the reasoning of the Appendix of [10]. Using the finite increments 
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Theorem, we have 

4>(tL)F'(f) - F'Q>(tL)f)\ < \Q>(tL) - I)F'(f)\ + \Q>(tL) - I)f 
So using similar arguments as previously, we get (with h = F'(f) and h = f) 



\III\ 



W2s-d/p,p 



< 



< 



1 ' (f)(tL)F'(f) - F'WL)f)) 2 tL4>(tL)f 



o 



9 J+\ V 2 



LP 



(4>(tL) - I)h tLcj)(tL)f 

1/2 



,-2s+d/p^\ 
t 



1/2 



LP 




where we used s < 2, which yields 
(14) 



tLcf>(tL)f 



< 



(tL) l ~ s l 2 (j)(tL)(tL) s l 2 f 



< t -d/(2p) 



(tLy/ 2 / 



< t -d/(2p)+s/2 



LP 



W 3 'P 



and the boundedness of the square functional associated to the function 



which is holo- 



morphic and vanishing at and at oo (see Proposition 1.13). We conclude the estimate of III 
since h = f or h = F'(f) belongs to W s ' p . The term IV is similarly estimated. □ 

Corollary 5.2. The diagonal term Rest (defined in Corollary 4-3) is bounded from LP x L q 
into U as soon as p,q G (1, oo] with < ~ := ^ + ^. Moreover for p € (1, oo), e > as small 
as we want and s > d/p (with s < (N — 2)/2), then 



\Rest(f,< 



\y/2s-d/ v ,p ^ 



iys+<!.p||5llW s + e >I>- 



Proof. Apply Theorem 5.1 to the quantities f + g and f — g with F(u) := u 2 . Then the 
polarization formulas give that Rest(/, g) has the same regularity has w in Theorem 5.1. □ 

Remark 5.3. Usually, we have a gain of regularity of order s — d/p for this quantity. Here we 
have a gain of s — d/p — e for every e > 0, as smal as we want. 

Remark 5.4. Of course, the assumption s G (d/p, 1) can be seen as very constraining. We 
want to explain here how it seems to us possible to weaken that point. First we point that this 
"technical" difficulty is new since it does not appear in the Euclidean situation. 

Legitimating the definition of the paraproducts (just before Definition 4-1 and Remark 4-%), we 
have developed (tL) 1 (<f>{tL) f '<j)(tL)F' '(f)) using the Leibniz rule of the Laplacian. Now for M « 
N, it is possible to do the same operation and develop (tL) M (4>(tL) f '4>(tL)F' '(f)) ■ We also obtain 
a sum of different terms involving higher order differential operators than previously. Similarly, 
in the beginning of the previous proof ( obtaining the several terms I, II, III, IV and V ), we 
have developed (tL) 1 F(cj)(tL)f). Taking the same exponent M, we can expand (tL) M F((j)(tL)f) 
in different quantities. Indeed, we find different multilinear differential operators Tjj such that 

(tL) M F((f>(tL)f) = T h M(tL)f, F'Q>(tL)f), F^" 1 ) (HtL)f)) 



+ F^ 2M \ ( t>(tL)f)\t 1/2 X4 ) (tL)f\ 



2M 
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By this way, we may define a new kind of paralinearization and prove that for f E W s+e ' p 
with large s > d/p then 

roc ^ _ _ jj. 

F (f) = E / T J , t (m)f,m)F\f),..-,m)F (2M - 1 Hf))-+w 

with w E W 2s - d / p ' p . By this way, the rest w should be decomposed as previously, in the corre- 
sponding term II, it will appear It^Xt^iL)/!^ -1 such that the corresponding square functional 
will be bounded as soon as {s — d/p){M — 1) > s — 1. The other terms may be bounded (as we did 
in Step 4) since they will have quantities as \(j){tL)F^ k \f) — F{k){4>{tL)f)\ and other differential 
operators on <j){tL)f. The key idea is that now the multilinearity of the operator Tjj will be 
sufficiently high to involve sufficiently such differential terms, each of them bringing a positive 
power oft as shown in (14)- 

By this way, it is also possible to get a paralinearization result for high regularity s > d/p 
and s « N (by taking a large exponent N), by defining new multilinear operators involving the 
derivatives F^ k \f). 

As explained in [10] (see its Appendix 1.3, theorem 38), a vector- valued version of the preceding 
result allows us to prove the following one: 

Theorem 5.5. Consider s E {d/p, 1), / E W s+k,p and a smooth function F{x, u\, ■ ■ ■ ujv) £ 
C°°{M x R N ) with F{x,0,--- ,0) = 0. Then by identifying {1, • • • ,N} with a set of multi- 
indices {ai, • • • «Ar} (and \ai\ <k), we can build 

(15) xeM ^F(x,X a J(x),---X aN f(x)) 
which belongs to W s ' p . Moreover, 

N 

(16) F(x,X a J(x),--- ,X aN f(x)) =^2H-[d Ui F](x,x ai f(x),- ,x aN f(x))(X a J)(x) + w{x) 

1=1 

with w E W 2s ~ d / p ' p . 

6. Propagation of low regularity for solutions of nonlinear PDEs 

As in the Euclidean case, paralinearization is a powerful tool to study nonlinear PDEs and 
to prove the propagation of regularity for solutions of such PDEs. Let us try to present some 
results in this direction with this new setting of Riemannian manifold. 

Let us consider a specific case of nonlinear PDEs for simplifying the exposition : let F(x, ui, ■ ■ ■ u K+ i) E 
C°°(M xR K+1 ) be a smooth function with F(x,0, ••• ,0) = 0. Then by identifying {1, •• • 
with a set of multi- indices {0, 1, • • • k}, we deal with the function 

(17) F(f,Xf) := x & M ^ F{x, f{x),X l f{x), - ■ ■ X K f{x)) 

for some function /. That corresponds to the case N = K + l, k = 1 with a\ = and = X^\ 
for t = 2,---JV + l in (17). 

Theorem 6.1. Consider s E {d/p, 1), / E W s+l,p and a smooth function (as above) F{x, U\, • • • un) E 
C°°{M x M. N ) with F{x, 0, • • • , 0) = and assume that f is a solution of 

F{f,Xf){x) = 0. 

Consider the vector field 

K + l 

T{x) :=^2[d u% F]{x,f{x),X l f{x),---,X K f{x))X i . 

i=2 

Then, locally around each point xq E M in "the direction T", the solution f has a regularity 
W s+1+P for every p > such that 

p < min{l, s — d/p}. 
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In the sense that 

U(f) := ^2[d Ui F](x,f(x),X 1 f(x),--- ,X K f{x))L^/ 2 XiU) G V ■ 



K + l 



i=2 



Such results can be seen as a kind of directional "Implicit function theorem", where the 
regularity of F(f,Xf) implies some directional regularity for / (in the suitable direction, where 
we can regularly "invert" the nonlinear equation). 

Proof. The previous paralinearization result yields that 

i=l 

which gives 

K 

T F (f) := £n^ +iW/) (L Q ^/) € L p , 
i=i 

where II is another paraproduct. Indeed 
U b (a) = - 



(tL) a i){tL) (tL) 1 - a (f>(tL)a<l>(tL)b 
{tL) a 4>(tL) \r a 4){tL)a4>{tL)b 



dt 

T 



dt 

T' 



where we have taken the notations of the definition for the initial paraproduct II (see Definition 
4.1). Then, we want to compare this quantity to the main one : U(f). So let us examine the 
difference. Since for every constant c, we have 



it comes 



hence 



cf = n c (/) = L a u c (L- a f) = n c (/), 

K 

u{f)(x) = ^2n [du ^ iF]{f{x))Xf{x)) (L a x l f)(x), 



i=l 



T F (f)(x) - U(f)(x) = ^U Kx (X l L a m 



i=i 



with \i >x (') = [d Ui+1 F](f, X f) — [d Ui+1 F](f(x),Xf(x)). It remains us to check that for each 
integer i, the function x — > Jl\. x (Xi(l + L) a f)(x) belongs to LP. Let us recall that 

dt 



U Xi jX t (l + Lrf)(x) 



(tL) a ^(tL) ^tL) 1 - a 4>(tL)X i L a f^>(tL)X i ^ (x)- 

r i dt 

(tL) a cf>(tL) [t- tt 4(tL)X i L a f<l>(tL)* i)ai \ (x)j. 



t 



Let us study only the first term / (the second one beeing similar): 

dt 



I :-- 



(tL) a ^(tL) (tL) l ~ a (f>(tL)XiL a f (f>(tL)\ i)X {x) 



< 



< 



1 



o J Mt x(B(x,t-V 2 )) 

POO f 

JO JCfx^H- 1 



1 + 



d(x,y) 



-d-8 



(tLy- a <l>(tL)X i L a f(y) ct>{tL)\ x (y) 



/ 2 ) 



(tLf-^t^X^fiy) 4>{tL)X i>x (y) 



dfi(y)dt 



dfi(y)dt 
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where we set C{x,2h~ x l 2 ) = B(x,T +1 t~ 1 / 2 )\B(x,2H~ 1 / 2 ) and by convention \C(x, 2H~ 1 / 2 )\ 
\C{x,2H^/ 2 )\. Now, for y € B(x, 2H- 1 / 2 ), we have 



1 + 



-d-<5 



y^ 2 jd-sk_f 



<E 2jd ~ Sh l 
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[d Ui+1 F]tf,Xf)(z) - [d Ui+1 F](f,Xf)(x)\ dfi{z) 



with i7 := [d Ui+1 F](f, X f) and C another systems of coronas. So we get 
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Using Cauchy- Schwartz inequality, we also have 
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where we have used a change of variables and 8 > d + 1. So using exponents q,r > p (later 



± = ± + i. 

other hand Fefferman-Stein inequality for the maximal operator, it comes 



chosen) such that - = - + ^ , boundedness of the square functional on the one hand and on the 
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Then, using the characterization of Sobolev norms (using this functional, see Proposition 3.6), 
we conclude to 

\\I\\lp ^ ||/||w 2 =>?||- H 1lw 1 > r - 
We also chose exponents q, r such that 

1 2a 1 8 + 1 Ills 
> — and - - - > - - -, 

q a p a r a p a 

which is possible since - < because of the condition on p. Then Sobolev embedding 

(Proposition 3.4) yields that W s+1 ' p W 2a > q and W s ' p <— > W 1 ^. Finally, the proof is also 
concluded since we obtain 

\\I\\lp < ||/||^+i,p||i?||w s 'P) 

which is bounded by / € W s+1 ' p (due to H := [d Ui+1 F](f, Xf) with Proposition 3.9). □ 

We let the reader to write the analog results for higher order nonlinear PDEs. 

Remark 6.2. Let us suppose that the geometry of the manifold allows us to use the following 
property: For a > 0, the commutators [Xi, (1 + L) a ] is an operator of order 2a, which means 
that for all p £ (1, oo) and s > 0, [X h (1 + L) a ] is bounded from W s+2a > p to W s ' p . 

This property holds as soon as we can define a suitable pseudo-differential calculus with sym- 
bolic rules : in particular, this is the case of H-type Lie groups, using a notion of Fourier 
transforms based on irreductible representations, see [5, 21]. 

Under this property, we can commute the vector field X with any power of the Laplacian and 
so with the same statement than in the previous theorem, we obtain that 

T(l-L)^ E fGL p . 

This new formulation better describes the fact that f S w s+p+1 ' p along the vector field V. 
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